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ABSTRACT 
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loss  is  quadratic.  In  this  paper  we  have  derived  certain 
optimal  properties  of  Stein’s  estimator  for  a  more  general 
loss  function.  It  is  shown  that  the  estimator  is  minimax 
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1.  Introduction .  Let  X  =  (X^,  X^)  be  a  vector  of  k 

independent  observations  from  k  normal  populations,  where 
X.  ^  N(e. ,1) ,  i  =  1,  ...,  k.  We  consider  the  problem  of 
estimating  the  vector  9  =  (9^,  9^)  .  Let 

6(X)  =  (5 L (X) ,  6k(X)) 

where  6^(X)  is  an  estimate  of  9^.  Let  the  loss  be  squared 
error,  given  by 

(1.1)  LC9.6)  *  ~  |  |«-0|  |2  =  i  l"  (5.-9  )2. 

K  K  i=l  1  1 

The  maximum  likelihood  estimator  5°(X)  =  (X^,  ....  X^)  has 
risk  equal  to 

R(9,6°)  -  Eq  LCS  ,6°)  =  1 

for  each  9.  James  and  Stein  (1961)  showed  that  the  estimator 
* 

5  (X) ,  given  by 

(1.2)  6*(X)  =  (l-ili)X. 

where  S  =  j  j  X !  j  w  ,  lias  risk  R(9,$  )<1  for  all  values  of  6, 

*  0 

if  k  _>  3 .  That  is,  6  (X)  dominates  5  CX).  Therefore,  the 

maximum  likelihood  estimator  is  inadmissible.  Even  though 

onlv  x.  seems  to  be  relevant  to  9.,  it  is  interesting  to  note 
l  i  ® 

* 

that  the  estimate  5^(X)  depends  on  the  entire  vector  of 
observations  X. 

The  surprising  result  of  James  and  Stein  has  stimulated 
considerable  research  in  the  last  two  decades,  on  the  problem 
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o£  estimating  the  mean  of  a  multivariate  normal  distribution. 
For  reference  to  some  recent  work  in  this  area,  see  Alam  (1975) 
and  Berger  and  Bock  (1976).  Almost  all  the  known  results 
assume  a  quadratic  loss,  chiefly  for  the  reason  that  nonqua¬ 
dratic  loss  is  difficult  to  handle  mathematically.  In  this 
paper  we  have  deviated  from  the  general  trend.  We  consider  a 
more  general  loss  function,  given  by 

Cl-3)  L  ( Q ,  6)  =  l  ^  (<5.  -  6.  )  2m 

K  i=l  11 

J> 

where  m  < is  a  positive  integer,  and  obtain  the  risk  of 
* 

<5  (X) ,  using  the  generalized  loss  function.  In  the  sequel  we 
shall  denote  by  R(0,o)  the  risk  of  an  estimator  6(X)  with 
respect  to  the  loss  (1.3). 

The  maximum  likelihood  estimator  has  a  constant  risk, 
given  by 

(1.4)  R(9,<SU)  =  T  (m+i)  . 

^  TT 

i  7 

On  the  other  hand,  R(9,6  )  is  a  function  of  ||9j|“  for  m  =  1. 

*  *  .  * 

For  m>l,  let  R  (A,o  )  denote  the  average  value  of  R(9,5  ) 

with  respect  to  the  distribution  of  9,  uniform  on  the  sphere 

(1.5)  ! 1 9 ! ! 2  =  X 

with  center  at  the  origin  and  radius  equal  to  A  .  In  the 
following  section  it  is  shown  that 

(1.6)  R*(a,5*)  .  R( 9 , 5°) 

for  all  values  of  \ ,  if  k  is  sufficiently  large. 
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The  proof  of  (1.6)  involves  very  difficult  computation 

It  is  perhaps  the  strongest  result  we  can  show  for  Stein's 

* 

estimator,  since  5  is  designed  for  a  spherically  symmetric 

loss,  whereas  (1.3)  is  spherically  non- symmetric .  On  the 

other  hand,  if  there  exists  an  estimator  which  has  smaller 
* 

risk  than  6  with  respect  to  (1.3),  its  risk  function  is 
likely  to  be  mathematically  intractable.  Therefore,  the 
inequality  (1.6),  which  is  the  main  result  of  this  paper, 
is  significant,  even  though  we  should  look  for  a  sphericall 
non- symmetric  estimator. 


i 

i 


■  >&}  J 1 


Let  m  =  1.  Consider  a  Bayesian  approach.  Suppose  that 


the  6^’s  are  independently  normally  distributed  with  mean 
zero  and  variance  A.  A  Bayes  estimator  of  9  with  respect  to 
the  given  prior  is  5(X)  where 

(1.7)  6 . I X)  =  (l-B)X.  ,  i-1,  ...»  k 

B  =  1/ ( 1 ♦A) . 

The  Bayes  risk  is  equal  to  1-B. 

If  B  is  not  known  then  S  is  .a  sufficient  statistic  for 

2 

B,  and  BS  is  marginally  distributed  as  >  chi-squared  with 
k  degrees  of  freedom.  Let  B(S)  be  an  estimator  of  B. 
Substituting  B(S)  for  B  in  (1.7)  we  get  an  empirical  Bayes 
estimator  r(X),  given  by 

Cl  .  S  )  •  .  t  X)  =  i  1-B  I'S)  )X  ,  i  =  1  ,  .  .  .  ,  k  . 

Let  .;  ( B ,  •  )  denote  the  Bayes  risk  of  <5.  The  difference  between 
(B,i  )  and  the  Bayes  risk  is  equal  to  B,  and 

,  ,  o  ,  -  1  B  >  ‘j  ,  -  ,  BjTbl)  2 

1  '  '  1  B  “ 1  B  ’ 

\ 

uhc  re  B T  "  ,  . 

k  -  1  ~  ■* 

Lor  B  l  S !  =  we  have  ;  i  X)  =  :  (XI,  and  the  left  side  of 

l..*)  is  equal  to  f.  Lfron  and  Morris  i  1 L1  - 3 )  have  shown  that 

K 
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the  minimax  value  of  the  right  side  of  (1.9)  is  equal  to 
2 

■£.  Thus  Stein's  estimator  is  minimax  in  an  empirical  Bayes 
sense,  considering  the  problem  of  estimating  B  with  loss 
function,  given  by 

L  (B ,  B)  =  (^)2  . 

A  similar  result  on  the  minimax  property  of  Stein's  estimator 
is  obtained  for  the  generalized  loss  (1.3). 

2. Main  results.  First  we  compute  a  Bayes  risk  of  the  esti¬ 
mator  6  (X)  given  by  (1.8),  using  the  loss  function  (1.3).  We 

make  the  Bayesian  assumption  ^  ~N(0,A)  as  a  convenient  mathe- 

•  d  d  7 

matical  tool.  Then  marginally,  X.  ~.N(0,  1+A)  and  BS  -  Xj)  . 

Let  r  be  a  non-negative  integer,  and  let 


(2.1) 


? 1  1 

d  -  4  T(r4) 

r  /IT 


(2.2) 


c  =  ~  E*(Q  Xjr)/Sr) 
r  K  i=  1  1 


1  *  4  7r  *  r 

=  r  (E  (E  XT  )/ (E  S r ) ) 

K  i=l  1 


2'r~  lb) 


rc4*n  r 


where  E  denotes  expectation  with  respect  to  the  marginal 


distribution  of  the  X.'s. 

i 


The  second  line  in  :2.2)  follows 


r*  ?  r  T 

from  the  fact  that  the  distribution  of  Q  X.  )/S  does 

i=  1  1 

not  depend  on  A,  that  S  is  a  complete  and  sufficient  statistic 

k  2r 

for  A,  and  hence  S  and  (2  X.  )/Sr  are  independently  distri- 

i=l  1 

buted.  Note  that  E  X2r  =  d  . 

i  r 

Conditionally,  given  X. 


(2.3) 


0i|Xi~  N((1-B)X.,1  -  B) 


Using  (2.3)  we  obtain  after  simplificat ion ,  the  Bayes  risk  of 
6 (X) ,  given  by 

-  m 


(2.4)  P  (B  ,5 )  =  l  c  d  ("‘i1)  (l-B)rE*(Sm_r(B(S)-B)2m'2r) 

r=0  m'r  r  ir 


m 


’m  2m*rr (T+m-n) 

I  Cm-r  dr  O  - IT -  B  (l'B): 

r=0  m  r  r  .r  r 


1 B(xk+2m-2r/B) 

i  B 


2m-  2r 


-  1 


-m  m-r  r  B  ^  xk+  7m*  ”*  B^ 

l  a  B  r  ( 1  -  B )  r  E| - k^-m-  ~r - 1 


r=0 


B 


_  ^  1  .  rn  n  ^  rn  ~  i*  k  \  y 

ar  =  cm-r  dr  !:r)  -  •  'T+tt-rJ/rt?) 


where 


V  (2r}  r(r+i}  rtra-r+I^ 


(m)  d  . 
r  m 


The  Bayes  risk  of  Stein's  estimator  is  obtained  by 
putting  B(S)  =  in  (2.4).  Then 


2m-  2r 


(2.5)  o (B,  5*)  =  l  a  Bm‘r(l-B)r 

r=0  X,  . 


k  +  2m-  2r 


Khere 


I  a  Bm*r(l-B)r  ( -  2m+  2r ,  -m+r;-  — ^ 

»=n  r  ~  - 


=  I  a  b  Bm~r(l-B) 


■  1  *  S *  *  £  * 


denotes  the  confluent  hypergeometric  function,  and 


(2.6) 


b  =  o(-2m+2r,  m+r;  -^-=-)  . 


Remark  1:  We  have  that  b  =1  and  b  -  0  as  k  -*  *  for  r<  m. 
-  m  r 

The  confluent  hypergeometric  function  J(a,b;x)  satisfies  the 
equat ion 


TnjI 


X 


d?  +  (b'x)  fl  -  ay  =  0 

The  above  equation  shows  that  for  a<0,  y>0  and  b  <_  x  <  0 ,  y 

is  a  convex  function  of  x  when  4^-  >  0.  Therefore,  if  v  =  0 

dx-  ’ 

for  x  =  -«  then  4^-  <0  for'all  x<0.  Hence,  b  is  decreasing 
dx  —  r  13 

in  k  for  r  <m,  and  therefore  there  is  a  minimum  value  of 
k  =  k^,  say,  depending  on  m,  such  that  b^_<l  for  r  =  0,1,..., m 
and  all  k>k  .  Numerical  computation  shows  that  k  is  the  mini- 
mum  value  of  k  for  which  bQ<_l.  Since 


C(B,6U)  -  l  a  Bmr(l-B)r  =  d  . 

r=0  T  m 


we  have  from  (2.5)  that 


(2.7) 


P  (B  ,  <5  *)  <  p  ( B  ,  5  ° ) 


for  k  >  k  and  all  values  of  B. 

—  m 

Let  V  4  X£(*0,  a  non-central  chi-squared  with  k  degrees 
of  freedom  and  non-centrality  parameter  X  =  :  !  6 j  ;  " .  The 
following  lemma  gives  a  formula  for  R  (.;,d  ),  the  average 
value  of  R  f  5 , 6  )  on  the  sphere  !  I  B  1  I  “  =  X. 


Lemma  2.1. 


(2.8)  R  ( ■'*  ,  <5  " )  = 


.k. 

rt-r)  m  ,  ..  .. 

-  l  (-l)ra  b  £(£(-!-, 4-m;4/(V)'m) 

r(T-n)  r=0  r  r  "  - 


m  r  . 

-  I  l-l)r3b  :  i - 1) r (^)^(m-t ,4; -4) 

r*0  t=0  L  ~  “  - 
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Proof:  Let  U  -x^/E-  Apriori,  9^  ?.N(0,A).  Therefore 

|[9|j"  ~  A  •  The  right  side  of  (2.8)  is  a  function  of 

|  | 9 |  | “  equal  to  g ( |  !  9  |  (  ^ ) ,  say.  Taking  expectation  with 

? 

respect  to  the  distribution  of  i|6||",  we  get 


(2.9) 

E  g  C  |  i  ©  !  ! 2 ) 


^2-*  rm  ,  nr  ,  -crrTr  k  UWU,-ITU 

— r -  l  (-1)  a  b  E((b(-r,--m;  j) (-)  ) 

r (|-m)  r=0  r  r  “ 

f  ab  Bm-r(l-B)r 
r=0  r  r 


=  p (B  ,  6  )  from  ( 2 . S)  . 

Since  e.  ~N(0,A),  the  conditional  distribution  of 
6  given  | { 9  (  [“  =  A  is  uniform  on  the  sphere  j  ' 0 J  |  ”  =  A 
There  fore 

(2.10)  t(B,6*)  =  E  R*(i  !  9  ||  ”  »<S*) 


where  the  expectation  on  the  right  side  of  (.2.10)  is  with 

*> 

respect  to  the  distribution  of  |  '  9  j  j  “  .  Since  the  distribution 
■) 

of  !  ]  9 ;  j  “  is  complete  for  B,  from  (2.9)  and  |2.1'1)  we  have 

(2.11)  R*(! i 9 ! !”,?*)  =  gC ; !e; i2) 

for  all  values  of  ! [ 9 |[ .  From  12.11)  we  get  the  first 
equality  in  (2.S).  The  second  line  in  (2.3)  is  obtained  from 
the  first  line,  using  the  following  formula. 

.  F  ifc-m+t  +  r) 

\  X  I 


\  /  -> 


-  10- 


.-X/2  r  *Vm+t)  -,k  k.  A, 

'  b-  0  l  7  y  7  9  7  3 

r(f)  “  z  ‘  z 


rtrm+t)  t,  „  k.  x, 

- p -  <P(m-t,  — ,  -tJ 

r  Cf)  22 


The  lemma  is  proved. 


Remark  2:  The  formula  (2.8 j  can  be  written  as 


Q.E.D. 


(2.13)  k 

*  *  _  ^  r  ^  tr^  ^ ~2 ^  V  -m+t 

R  (A, 5  )  =  l  (-l)ra  b  [  C-D  O  — l - 

r=0  r  r  t=0  t  r(|-m+t)  Z 


where 


(2.14) 


„ ,V. -m+t 
et 


,  ,st+t  fr,  1 

1  (-D  arbr  ^t  — k - 

r=t  r  r  1  r(-|-m+t) 


values  of  et  for  t  =  m  and  m  -  1  are  computed  easily. 


he  have  e  =  d  and 
m  m 


(k- 2) “ 

—  n.  -  mdm 
2k  m 


Note  that  the  risk  of  the  maximum  likelihood  estimator  is 
equal  to  dm  and  that  for  m  =  1,  that  is  for  squared  error  loss, 


we  have 
(2.15) 


R(e,6*)  =  R* ( |  !  6  !  !  2 ,5*) 


,  (k  -  2)  "  F  1 
~  k  L  V 


0 

-  Px(e,5  )  . 


m 


T 


-**■ 
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We  shall  now  obtain  the  main  result.  Let  m>l. 

From  the  integral  representation  of  the  confluent  hyper¬ 
geometric  function,  given  by 


±(a,b, ;x) 


=  r(b) 

r (a)  r (b-a) 


1 

J  0 


xu  a-1,,  .b-a-1, 

e  u  (1-u)  du 


b>  a>  0 


—  k  a 

we  find  that  the  value  of  ^_(m- 1  lies  between  0  and  1, 

and  that  it  tends  to  1  as  k  ■+■  °° .  Since  b  =1  and  b  +  0  as 

m  r 

k-*-oo  for  r  <  m,  from  (2.8)  we  have  that 

is  is 

(2.16)  R  (\,5  )  -*•  0  as  k  , 


On  the  other  hand,  as  A  the  summation  in  the 
on  the  right  side  of  (2.13)  is  dominated  by  the 
terms,  equal  to 


(2.17) 


e  + 
m 


e 

m 


1 


-1 


(k-2)2 

k 


second  line 
sum  of  two 


Hence 


<  d  (1- 
m 


m (k- 2 )  " 
k  (k+  ) 


)  • 


(^•18)  R*(A,6*)  dm 

for  sufficiently  large  values  of  A. 

Next  we  consider  the  case  in  which  both  k  and  A  are  large. 

Expanding  b  ,  given  by  (2.6),  in  negative  powers  of  k  we  find 

->  -  2 

that  b  =  1,  b  .  =  -p  and  b  =  0(k  )  for  r<m-l.  Therefore, 

m  m-1  k  r 

for  large  k,  the  summation  on  the  right  side  of  (2.8)  is 
dominated  by  the  sum  of  two  terms,  equal  to  D,  say  where 


(2.19)  D 


m  r  -r  r  r  -  k  A 

V  i  -  r  a  b  :■  (.  - 1 ) 1  (5)Mm- 1 -4) 

“r=m- 1  1  r  “t-0  "  “ 
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-X/2 


m 


e  l  C-Dra  b  f 

r=m-l  t=0  c  “ 


r  (f  mn+T)  P(|) 


m  r  i  (.y-m+t+i; 

El  (-Dra  b  l  C-l)tCj)  — | - r 

r=m- !  r  r  t-o  r  r(4+T)  T(| 


m+t ) 


m 


■  E  Iri>.HlVr  Jt,0('1)t(t)(Trl77)m't(1'rra7TTTtct 


+  0(k'^)) 

where  T  is  a  random  variable  distributed  according  to  the 
Poisson  distribution  with  mean  equal  to  A/2  and  ct  = 

(m-t) (m-t+l)/2.  After  simplification  we  get 


m-  It-  — 

j  T+2m  m(m+l)k/2 

=  d  E 
n 

T 

T+k/2 

l  J 

j^T+k/2  2  (T  +  k/2)  2  _ 

-  A  p  1 

/  > 

T 

m'1gk(T)  +  0(k'2) 

-  Cl  C 

m 

T+k/2 

+  0(k'2). 


where  g,  (T)  <  1  for  k  >  4.m.  Now 

K 


T+k/2 


m- 1  <  £  |  T 


—  j  T+k/  2 

k  J 

k/2 


1  1  - 

=  1  - 
<  1  - 


ET+k/2 

k 


A+  k 

1 

l+k 


for  .\  <  k * 


*  * 

R  )  <d 


m 


There  fore 
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for  sufficiently  large  k  and  X  _<  k  .  On  the  other  hand,  from 

(2.13)  and  (2.17)  it  follows  that  the  above  inequality  holds 

2 

also  for  fficiently  large  k  and  X  _>  k  .  Therefore,  we  have 
the  following  theorem. 

2 

Theorem  2.1.  Let  X  =  |  j  © |  j  .  Given  a  positive  integer  m, 

* 

there  exists  a  value  of  k  =  k  ,  say,  depending  on  m,  such 

that,  R  ( X ,  6  )  <d  =  R(0,6^)  for  all  6  and  k  _>  k  . 

m  m 

Next  we  show  that  Stein's  estimator  is  minimax  in  an 


empirical  Bayes  sense.  From  (2.4)  we  have 


(2.20)  p(R,6)  =  d  l  (™) Br ( 1- B) 1 

m  r=0  r 


B  (W  ) 


where  W  ~  x£+Tr/B.  In  deriving  the  above  result  we  have 
assumed  that  9^  ~N(0,A).  We  can  stop  short  of  full  Bayesian- 

hood  by  assuming  that  A  or  equivalently  B,  is  unknown  and  must 
be  estimated.  That  is,  we  become  empirical  Bayesian,  considering 
B(Wf)  as  an  estimate  of  B. 

Consider  a  beta  prior  distribution  on  B  with  density 
function 


g(B)  =  3s" 1 (1- B) 


n-l  r rs+n) 

ru)r(n) 


4  >  0  ,  n  >  0 . 


fl(B,6)  =  (D  (  B  ,  <5 )  -  ( 1  -  B)  a  1  /B 

m 


Let 
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From  (2.20)  we  have 

(2.21) 

rl  ~  m  rl  m  m  f  B  (W  ) 

I  Q(B.  <S)g(B)dB  >  l  d  I  (™I) Br~  i(l-B)m  r  — V- 

J o  J'o  r  <•  B 


1  g(B) dB 


where 


L  B2  "  (l-B)m'r+n'1e'BWr/2dB 
B(Wr)  =  - 

f0  B“  ’  (l-B)m'r+n*1e'B'Vr/2dB 


•e  r'  “dB 


(|+M)_K^£;,T+m-r*c+n;  ~V2) 

(y+m -  r  +  4  +  n -  1 )  0  (j*  £,  - 1 ,  y+m -  r  +  £+  n '  1 ;  -  wr /  2 ) 


From  the  asymptotic  property  of  the  confluent  hypergeometric 
function  (see  e.g.  Erdelyi  (1955)  §6.13.1)  it  is  seen  that 


N  B  (w_)  k+,c.2 

(2.22)  v  ( W r )  =  -^-k-,-v-r  -  K  ■“-?/  as  Wy  -  ». 


Let  %-0.  Then  B  converges  in  distribution  to  :ero  and 
Wr  converges  in  distribution  to  *>.  Therefore,  v  (W  )  converges 
in  distribution  to  1,  and  from  (2.21)  we  have 
(2.23) 


Sup.'.fB,  5)  >  Lim  I  T(B,ng(B)dB 

R  r-0  1  0 


■1  m  r-1  ,  _ 

rT  D>m-rr,k-:„ 


I”  dmj0^r=1VB  (l-B)m-rE(|^(Wr)-l)-r  g  (B1  dB 
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/ 


1 1  m 
= Lim  d  | 


(’!)Br(l1-B)m'rE(|f  -1J 
0  mj0  r-1  r  B  r 


g(B)dB 


I1 

= Lin  [  IHB,  5  jg (B) dB 

s-o  J  0 


a  .  b  . 
m  -  1  m- 1 


*  *  * 

=  Sup  "2(3,6  )  for  k  2.  k  ,  say. 
B  m 


The  last  equality  in  (2-23)  follows  from  the  fact  that 
_  1 

bf  *  0(k  “)  for  r<m  -  1.  Thus  we  have  obtained  the  following 
result  from  the  empirical  Bayes  approach,  considered  above. 

Theorem  2.2.  The  minimax  value  of  the  Bayes  risk  o(B,5)  is  * 

■>m  *  ** 

equal  to  J  =  Sup  f.  ( B  ,  5  )  for  k  >  k 

K  O  g  ITl 

Remark  5 :  Clearly,  the  results  given  above  are  appli¬ 
cable  to  the  case  in  which  the  loss  is  a  polynomial  function 
of  the  squared  error  with  non-negative  coe f f ic ient s .  It 
might  be  possible  to  extend  the  given  results  to  the  case  in 
which  the  loss  is  a  monotone  increasing  convex  function  of 
the  squared  error. 

it  *  * 

Table  I  below  gives  values  of  k  and  k  for  m  =  1(1)10. 

-  m  m 

* 

For  the  given  values  of  m  it  is  found  that  k  =  k  where  k  is 

m  m  m 

the  smallest  value  of  k  for  which  (!.")  holds.  Table  II  below 

*  *  A 

gives  values  of  R  ('•.,5  )d  for  n  =  1 1  1 )  1 0  and  *.  and  k  *  k  . 
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TABLE  II  -  Values  of  R*(A,6*)/d  for  k  *  k* 

^  ’  1  m  m 


ss 

B 

2 

3 

B 

5 

6 

B 

8 

9 

10 

0 

.  6666 

.  5714 

.  8581 

.  5301 

.  8640 

.  5758 

.9178 

.  6247 

.  977  9 

.6"24 

1 

.7534 

.  5294 

.6444 

.4004 

.  6176 

.4202 

.6532 

.4522 

.6959 

.4851 

2 

.8206 

.  5239 

.5137 

.3140 

.4513 

.3119 

.4709 

.  3311 

.  5000 

.3531 

5 

.  8634 

.  5379 

.4377 

.2573 

.  3385 

.2360 

.3444 

.2455 

.  3631 

.  2594 

4 

.  8933 

.  5615 

.  3976 

.2212 

.  26  20 

.  1826 

.  2560 

.  1345 

.  2666 

.1925 

5 

.9145 

.  5889 

.  3808 

.  1997 

.  2101 

.1447 

.  1937 

.  1407 

.1981 

.  1444 

6 

.  9299 

.6170 

.  3790 

.  1886 

.  1753 

.1179 

.  1494 

.1089 

.  1491 

.1095 

n 

.9412 

.  6440 

.  3866 

.  1849 

.1525 

.0989 

.1178 

.0858 

.1158 

.  0841 

10 

.  9614 

.7132 

.4357 

.  2008 

.1264 

.0703 

.  0669 

.  0468 

.0554 

.0410 

15 

.  9757 

.  7899 

.  5302 

.  2640 

.  1454 

.0672 

.  0442 

.  0257 

.  0236 

.  0165 

20 

.  9823 

.8361 

.  6066 

.3346 

.1899 

.  0863 

.  0477 

.0231 

.  0161 

.009"  | 

25 

.  9860 

.8661 

.  6642 

.  3995 

.2405 

.1152 

.  0618 

.  0278 

.0160 

.0082  ! 

1 

100 

.  9966 

.  9642 

.  3958 

.  7640 

.  6397 

.4744 

.  3515 

.  2207 

.1410 

. 0" 20  ! 

1 

1000 

1 

\ 

.  9996 

.9964 

.  9893 

.9746 

.  9589 

.  9340 

.  9105 

.8769 

.8468 

.8061  ' 

l 

1 
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